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Abstract 

The gradient-flow dynamics of an arbitrary geometric quantity is derived using a general- 
■ ization of Darcy's Law. We consider flows in both Lagrangian and Eulerian formulations. 

The Lagrangian formulation includes a dissipative modification of fluid mechanics. Eu- 
lerian equations for self-organization of scalars, 1-forms and 2-forms are shown to reduce 
to nonlocal characteristic equations. We identify singular solutions of these equations 
corresponding to collapsed (clumped) states and discuss their evolution. 
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1 Introduction 



Singular solutions in continuum models. Under certain conditions, advection equa- 
tions in continuum dynamics support singular 5-like solutions. The simplest example is the 
conservation equation for the mass density p(x, t) of the form 

d t p + div(pu) = . 

As one may verify, this equation admits the single-particle trajectory 

p(x,i) = io<5(x-Q(t)) 

as its solution, for a constant weight w, provided the velocity vector field satisfies u(Q(t), t) = 
Q(t) and is sufficiently smooth. 

More generally, singular solutions exist for any tensor quantity undergoing the flow of a 
smooth vector field. One writes 

d t T + £ U J = , 

where T is a generic tensor field and £ u stands for the Lie derivative with respect to the vector 
field u, so that the corresponding term represents fluid transport [TJ. The singular solutions 
for this equation may be written in general as 

N 

T(x, t) = J2 Pifa *) " Qi(*, 0) ^ , 

where the sum refers to iV different singular solutions and s is a coordinate on the 2-th 
submanifold in M. 3 (a curve or surface). A famous example of this more general case is 
provided by Euler's vorticity equation for ideal fluids: 

dtUJ = —£ u (jJ = curl(u; x u) , 

with nonlocal velocity u = curl" 1 uj. In this case the tensor quantity U3 is an exact two-form 
and the singular vortex solutions for this system are supported either on a curve (in the case 
of filaments) or on a surface (in the case of sheets) in R 3 [2]. In two dimensions, one has point 
vortices instead. 

Emergent singularities in hydrodynamics. Another interesting example of this kind 
of advection equation is provided by the dispersionless Camassa-Holm (CH) equation in one 
dimension [1] or its generalization in higher dimensions, the EPDiff equation [3]. The dis- 
persionless CH and EPDiff equations possess singular solutions that emerge spontaneously 
from any smooth initial distribution. In particular, these are advection equations for the fluid 
momentum represented in geometric terms by a one-form density m <8> dV = m ■ dx ® dV. 
Just as in the case of vorticity dynamics, the momentum conservation equation 

(d t + £ u \ (m ® dV) = 

involves a nonlocal velocity vector field u = G * m, where G is the Green's function for a 
positive definite symmetric operator, typically chosen to be G = (1 — A) -1 . The singular so- 
lutions of CH and EPDiff emerge spontaneously from any confined smooth initial distribution 
of velocity and their pairwise collision interaction is completely known. (In the case of the 
CH equation in ID these solutions are solitons [3].) 
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Emergent singularities in aggregation dynamics. In recent years, another class of 
dissipative continuum models has been found to possess emergent singular solutions [5J, E]. 
These equations appeared in the study of aggregation and self-assembly phenomena at different 
scales in physics and biology. The main idea is that the approach to a critical point in 
free energy of a continuum material may concentrate matter into self-organized structures 
(clumps). Diverse examples of such processes include the formation of stars, galaxies and 
solar systems at large scales, growth of colonies of organisms at mesoscales and self-assembly 
of proteins, nanotubes or micro/nanodevices at micro- and nanoscales [TJ. Some of these 
processes, such as nano-scale self-assembly of molecules, are of great technological interest. 

Mathematically, the classic examples of continuous equations for aggregation are those of 
Debye-Huckel [8] and Keller-Segel (KS) [9] for the density of particles. For a recent review 
of developments in this area with an emphasis on biophysical modeling, see, for example 
[TO] . When diffusion is negligible, the physics of these models consists of a conservation law: 
d t p + div (p u V$) = 0, in which the "mobility" a may also depend on the density, as it 
happens for the collective potential $[/?]. This simple relation is known as "Darcy's Law" 
[TTj . Previous investigations have extended Darcy Law to incorporate nonlocal, nonlinear and 
anisotropic effects in self-organization of aggregating particles of finite size [12] . 

Besides the spontaneous emergence of singular solutions, the equations for aggregation 
dynamics exhibit self-assembly through the interaction among their singularities after they 
have emerged. Indeed, the study in [5J E] has shown how the singular solutions clump 
together after their formation. This feature reflects the physics underlying the mathematical 
description and distinguishes this class of dissipative equations from the energy-conserving 
equations in fluid mechanics, such as Euler's equation or the EPDiff equation mentioned 
above. 

Geometric gradient-flow equations. The goal of this paper is to formulate a principle 
for deriving and analyzing evolution equations for continuum systems that may exhibit self- 
assembly under a flow which reflects the mathematical properties of Darcy's Law, although 
involving an appropriate physical quantity possessing other geometric properties than particle 
density. A useful concept for deriving a continuum description of macroscopic pattern for- 
mation (e.g., aggregation) due to microscopic processes is the notion of order parameter. 
Order parameters are continuum variables that describe macroscopic effects due to micro- 
scopic variations of the material structure. They take values in a vector space called the order 
parameter space that respects the underlying geometric structure of the microscopic variables. 
A familiar example is the description of the local directional asymmetries of nematic liquid 
crystal molecules by a spatially and temporally varying macroscopic continuum field of un- 
signed unit vectors called "directors, see, e.g., Chandrasekhar [13] and de Gennes and Prost 

In order to establish macroscopic continuum models for aggregation phenomena of order 
parameters, the present treatment formulates such equations by generalizing the geometric 
structure underlying Darcy's law. This is a non-trivial step, since Darcy's Law possesses an 
ambiguity when written in a covariant form using the Lie derivative, as explained in the first 
part of the paper. This ambiguity may be eliminated by requiring the existence of singular 
solutions in the geometric gradient-flow (GGF) equations for the order parameters. The 
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requirement of singular solutions is physically justified by the same arguments that hold for 
Darcy's Law, that is recovering self-assembly of many-particle systems in terms of coherent 
structures. Self-assembly will have a specific physical meaning, depending on the nature of 
the particular order parameter. For example, it will represent alignment-like phenomena in 
the case of rods or directors. 

Plan of the work. Section 12.11 presents the geometric structure of Darcy's Law and the 
related ambiguity in its underlying geometric structure. Section 12.21 identifies a unique ge- 
ometric class of order parameter equations by requiring that singular solutions exist, so as 
to capture coherent structures by introducing a spatial averaging that follows them in a La- 
grangian sense. In Section 12.31 this Lagrangian averaging approach yields nonlocal equations 
expressed in fixed spatial coordinates. Rather general mathematical and physical requirements 
(correct geometry of the measured variable and independence of the method of measurement) 
in Section 12.51 yield again the allowable equations discussed in Section 12.11 Remarkably, 
these general considerations yield the same evolution equation (|18[) for an arbitrary geometric 
quantity as found by using the Lagrangian averaging method. The ambiguity in the roles of 
the order parameter and its mobility is overcome by focusing on singular solutions resembling 
aggregation phenomena. 

In Section [3] we develop a form of dissipation in terms of Leibnitz-bracket structures for for 
fluid mechanics, which reflect the Darcy dissipation term. As we show, this formalism leads 
to an interesting equation for modified dissipative fluid dynamics that does not require extra 
boundary conditions (unlike viscosity dissipation in Navier-Stokes equations). The formalism 
generalizes earlier modified fluid equations of this type in Bloch et al. [151 EE] an d Vallis et 
al. HZ]. 

Next, in Section [4] we consider a self-organizing system of an arbitrary geometric quantity. 
The principle for obtaining such evolution equations applies generally for any geometric type of 
continuum physical quantity (density, scalar, 1-form, 2-form, etc.). Thus, we obtain a family 
of evolution equations for physical quantities such as active scalars, momenta and fluxes 
(Sec. [5]) that nonlinearly influence their own evolution. These quantities convect themselves 
by inducing a velocity appropriate for continuum motion of any geometric order parameter. 
In other words, they can be written as nonlocal characteristic equations, as shown in Section 

SI 

After the geometric gradient-flow (GGF) equation 018p has been derived from general 
principles, we discuss its most remarkable feature; namely, this equation possesses singular 
solutions. In these singular solutions, each type of order parameter may be localized into 
delta functions distributed along embedded subspaces moving through the ambient space. 
Depending on the geometric type of the order parameter, the space of singular solutions may 
either form an invariant manifold, or these solutions may emerge from smooth confined initial 
conditions. In the latter case, the singular solutions dominate the long-term aggregation 
dynamics. From the physical point of view, such localized, or quenched solutions would 
form the core of the processes of self-assembly and are therefore of great practical interest. 
The formation of these localized solutions is driven by a combination of nonlinearity and 
nonlocality. Their evolution admits a reduced description, expressed completely in terms of 
coordinates on their singular embedded subspaces. 
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2 Geometric gradient-flow equations 



2.1 Geometric structure of Darcy law 

The Darcy law for the geometric order parameter p (density) [5j [6] is written in terms of an 
energy functional E = E[p] and a mobility p which takes into account of the typical size of 
the particles in the system (in general it depends on p). In formulas, one has the equation 

a " divUMv^y (i) 



dt —y-^-sp 

This may be stated by using the Lie derivative operator in the two possible ways, 

. dp n r . dp n . . 

either — = Jb , j u p or — = Jb , p , (2) 

where sharp (•)* denotes raising the vector index from covariant to contravariant, so its 
divergence may be taken (the sign in the right hand side is taken in agreement with the 
dissipative nature of the dynamics, as it is shown in Sec. 12. 4(1 . The evident difference between 
these two forms is that, unlike the first form, the second equation can be written as the 
characteristic equation 

^(x(t),t)=0 along ^ = u [p}= L[p]V 5 -^) (3) 

so that velocity u depends on density p through the gradient of the variation of the free energy 
E (velocity proportional to thermodynamic force with mobility p[p\) [51 161 [12] . 

As mentioned in the introduction, we want to generalize this geometric flow method under- 
lying Darcy Law to apply to other order parameters (denoted by k) with different geometrical 
meaning (not just densities). The key question for understanding the physical modeling that 
would be needed in making such a generalization is, "What is the corresponding Darcy Law 
for an order parameter «?" Namely, how does one determine the corresponding geometric 
flow for an arbitrary geometrical quantity k? The first problem is that there is no reason to 
consider only one of the two geometric formulations in (T5]). Although a characteristic form 
would be preferable because of its richer geometric meaning, no choice can be performed a 
priori. 

As a further step in the investigation of the geometric structure in Darcy law ([T]), one 
wants to write the corresponding equation ([I]) in a fully covariant fashion. To this purpose, 
one defines the diamond operator o as the dual of the Lie derivative under integration by 
parts for any pair (k, b) of dual variables and any vector field v [12]. That is 

(ko6,v) = (K,-£ v b). (4) 

In this way the diamond operation is defined as o : V x V* — > £*, where V is the vector space 
which k belongs to and V* is its dual, while X* denotes the dual space of vector fields in M 3 . 
In the case of Darcy Law, k = p is a density variable and b is a scalar function, so that one 
has 

( /0O 6,v) = -(p,vV6) = -(pV6,v). 
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Thus, by using the definition of diamond, one expresses Darcy Law ([I]) in covariant form as 



dp dp 
either -^r + £ / xe ,n» lAp] = or ^7 + ^/ XB .\t P = Q • ( 5 ) 



'pof> 



This formulation is not only covariant, but it can as well be extended to any quantity for 
which a Lie derivative can be defined, i.e. any tensor field, or geometric order parameter 
k. It only suffices to substitute p with k in order to obtain the following two possibilities: 

d Kj dtz 
either — = —£, ap ,h u\k] or — = —£, . . xw *tK. (6) 

It is important to notice that, unlike Darcy 's Law, the equations above are not invariant with 
respect to the substitution k <-> //[re] and the two equations in (jB]) are different equations. 



2.2 GGF equations and singular solutions 

At this point it is necessary to make a definitive choice between these two possibilities, although 
both appear to be consistent geometric extensions of Darcy Law (JT]). Actually, one recalls 
that equation §T§ with generic mobility p = p[p] has one more feature, besides its purely 
geometric character. This feature is the spontaneous emergence of singular solutions, 
that correspond in one dimension to the trajectories of N interacting particles in the system 
[5j [6]. This point leads to the question: is it possible to generalize the existence of singular 
solutions to geometric gradient-flow (GGF) equations? Thus one is motivated to look at one 
of the forms in and see whether they have this property. For example, in one dimension 
one looks for solutions of the form 

N 

K{x,t) = J2Pi(t)S(x-Qi(t)). 

1=1 

Upon pairing the first equation in ([6]) with a dual element 0, direct substitution of the singular 
solution ansatz yields 

{% *) = E fr ■* «<<*» + E w ■ *' «'<«» 

= - { £ (, o t r ' ) = - ° 0> ( k ° f ) ) = - ( k ' ^ o ^« f ) 

= - y Pi (t) _i 4 °-f- (7) 

where the symbol _l denotes contraction of indexes. In order for the singular solutions to 
exist, one would match terms in <fi and <$' and obtain the evolution equations for p n and Q n , 
as it happens for the density variable p in Darcy law [3 [6]. However, in general the term in 
the last line may involve higher derivatives, not just first order (for example, if k is a one-form 
density, then diamond is again a Lie derivative, which generates second order derivatives in 
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cf)). Therefore, the first choice in §5§ is not suitable to recover the singular solutions in the 
general case of an order parameter k. Instead, by following the same procedure for the second 
equation in (J6]), one obtains 



x=Qi(t) 



Now, since the Lie derivative is a first-order differential operator that is linear in the deriva- 
tives, one recognizes that the last term on the second line contains only terms that are linear 
in and its first order derivatives, while it does not involve any of its higher order derivatives. 
Thus, in higher dimensions one finds the following conclusion 

Theorem 1 The second equation of (Tj|) always allows for singular solutions of the form 

N 

«(x, t) = J2 5 ( x - Q*( s > *)) ds (9) 

i=l J 

for any tensor field k, provided \i and 5E/5p are sufficiently smooth. 

The proof of this statement is a verification through direct substitution, similar to the calcu- 
lation in equation (jSJ). As in the introduction, the variable s is a coordinate on a submanifold 
of R 3 : if s is a one-dimensional coordinate, then k is supported on a curve (filament), if s is 
two dimensional, then k is supported on a surface (sheet) immersed in physical space. 

At this point, one defines geometric gradient-flow ( GGF) equations as characteristic 
equations of the type 

5 w«),o = o aio„ g Tt = {^ TK ) HO) 

or, in Eulerian coordinates, 

where the quantity (// o SE/Sk)^ can be called generalized Darcy velocity and in the case of 
Darcy Law expresses the fact that the macroscopic velocity is assumed proportional to the 
collective force v oc V$, where the coefficient of proportionality is given by the mobility a. 
The fact that Darcy law ([T]) is invariant with respect to permutations of p and p is the reason 
why singular solutions are recovered for both possibilities in (j3J). This property is peculiar of 
Darcy law and does not hold in general. The distinction between the two cases identifies the 
geometric structure of the GGF family of equations. 

Variational formulation. As a further step one seeks a variational formulation. One 
defines a gradient flow by setting the L 2 pairing of equation fill I) with a test variable <fi (dual 
to k) equal to the variation SE of the free energy 

Ok A /, 8E\ 

OK, — ) , 12 



dt ' / \ '5k 
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where the variation 5k satisfies 



£ 




it At[«] 



(13) 



in order to recover equation ffTTj) . Indeed, one obtains 



( 



Ok 
dt' 



> 




) 




(14) 



The next question in the formulation of GGF theory is the particular meaning assumed 
by the generalized mobility [i[k\. This quantity has been related to the typical particle size in 
Darcy dynamics [3 [6], but the physical meaning of this quantity is not yet clear in the case 
of a generic GGF equation for the order parameter k. The next section presents the mobility 
as a smoothed quantity that keeps into account the dynamics of jammed states in the system, 
by introducing a typical length-scale. 

2.3 Motivation for multi-scale variations 

We seek a variational principle for a continuum description of coherent structures. This 
includes the evolution of particles of finite size that may clump together under crowded condi- 
tions. In crowded conditions, finite-sized particles typically reach jammed states, sometimes 
called rafts, that may be locally locked together over a coherence length of several particle-size 
scales. Thus, a variational principle for the evolution of coherent structures such as jammed 
states in particle aggregation must accommodate more than one length scale. A multi-scale 
variational principle may be derived by considering the variations as being applied to rafts, 
or patches, of jammed states of a certain size (the coherence length). For example, one might 
introduce a probability distribution for the sizes formed in the process of particle clumping. 
We shall take a simpler approach based on applying a Lagrangian coordinate average that 
moves with the clumps of particles. In this approach, the variation (8k) of the order parameter 
(k) at a given fixed point in space is determined by a family of smooth maps </?(s) depending 
continuously on a parameter s and acting on the average value k defined by 



in the frame of motion of the jammed state. Thus, y is a Lagrangian coordinate, defined in 
that frame. This motion itself is to be determined by the variational principle. The kernel H 
represents the typical size and shape of the coherent structures, which in this example would 
be rafts of close-packed finite-size particles. The mobility fi of the rafts depends on k, and so 
the corresponding variation of the local quantity k at a fixed point in space may be modeled 
as 



Here (fi(s)y = x(s) is a point in space, which tp 1 (s) returns to its Lagrangian label y and (p{0) 
is the identity operation. The average k = H * k is applied in a Lagrangian sense, following a 




(15) 




(16) 
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locally locked raft of particles along a curve parameterized by time t in the family of smooth 
maps. The latter represents the motion of the raft as ip(t)y = x(t), whose velocity tangent 
vector is still to be determined. When composed from the right the derivative at the identity 
of the action of (p(s) results in a variation 5k at a fixed point in space given by 

6k = -£ w { v )h[k\ with v(<p) = <^V _1 | S=0 > ( 17 ) 

where £ v ^^[k] denotes the Lie derivative of the mobility [/,[k] with respect to the vector field 
v(ip). In what follows we will treat /j,[k] as a general functional of k, not just a function of R. 

Taking the steps for k corresponding to those for p in ffl4|) but with the Lie derivative 
appropriate for the geometric nature of k allows one to close the GGF equation ffTTT) explicitly 

as 

-^ = -£u{ K ]K, in which u[k] = (jjl[k] o— J . (18) 

Here the co- vector (//[k] o6E/6k) resulting from the diamond operation is defined as in equa- 
tion (jlj) and sharp ( ■ )" raises its vector index to make it contravariant. Although Lagrangian 
averaging was invoked in defining the mobility as a coherence property, equation (1181) is ex- 
pressed in fixed spatial coordinates. Its Lagrangian heritage is recognized upon rewriting it 
equivalently in the characteristic form ([3]). 



2.4 A dissipative bracket for GGF equations 

We shall use the method of multi-scale variations to derive several versions of the geometric 
gradient-flow (GGF) equation ( TTTT) for different geometrical types of physical quantities k. In 
each case, the result is a novel nonlocal characteristic equation which possesses localized (sin- 
gular) solutions. Before developing these results, however, we shall first address the question 
of energy dissipation. This will allow us to introduce the dissipative bracket and to write 
the equations in an alternative bracket form. The corresponding energy equation follows from 
(ED as 

dE idK 8E\ I 6E\ /( SE\ ( 8E\i\ 

The formula for energy in (1T91 suggests the following bracket notation for the time derivative 
of a functional F[k], 

dF[K] _ /8k 5F\ I 8F\ 

The properties of the GGF brackets {{ E , F }} defined in equation fl20l) are determined by the 
diamond operation and the choice of the mobility /a[k]. For physical applications, one should 
choose a mobility that satisfies strict dissipation of energy, i.e. {{ E , E }} < 0. A particular 
example of mobility that satisfies the energy dissipation requirement is = kM[k], where 
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M[k] > is a non-negative scalar functional of k. (That is, M[k] is a number.) Requiring 
the mobility to produce energy dissipation does not limit the mathematical properties of the 
GGF bracket. For example, the dissipative bracket possesses the Leibnitz property with any 
choice of mobility. That is, it satisfies the Leibnitz rule for the derivative of a product of 
functionals. In addition, the dissipative bracket formulation (120]) allows one to reformulate 
the GGF equation fill I) in terms of flow on a Riemannian manifold with a metric defined 
through the dissipation bracket, as discussed in more detail in |12j . 



Previous dissipative brackets. Historically, the use of symmetric brackets for introducing 
dissipation into Hamiltonian systems seems to have originated with works of Grmela |18j . 
Kaufman [19] and Morrison [20] . See [21] for references and further engineering developments. 
This approach introduces a sum of two brackets, one describing the Hamiltonian part of the 
motion and the other obtained by representing the dissipation with a symmetric bracket 
operation involving an entropy defined for that purpose. Being expressed in terms of the 
diamond operation ( o ) for an arbitrary geometric order parameter k, the dissipative bracket 
in equation (120]) differs from symmetric brackets proposed in earlier work. The geometric 
advection law (|T8l) for the order parameter will be shown below to arise from thermodynamic 
principles that naturally yield the dissipative bracket (120]) . Moreover, being written as a Lie 
derivative, the equation of motion (1181) respects the geometry of the transported quantity. The 
dissipative brackets from the earlier literature do not appear to be expressible as a geometric 
transport equation in Lie derivative form. 



2.5 Thermodynamic and geometric implications 



Equations (fT81) may be justified using general principles of thermodynamics and geometry. 
Consider using an arbitrary functional F in (120]) as a basis for the derivation of an equation 
for k. Suppose k is an observable quantity for a physical system, and that system evolves 
due to the inherent free energy E[k] in the absence of external forces. This is the physical 
picture one envisions, for example, when thinking about any process of aggregation. Suppose 
we would like to measure the time evolution of a functional F[k], which may for example 
represent total energy or total momentum. For an arbitrary functional F[k] and for a given 
free energy E[k] one finds, in general 



dF 



d t K 



5_F_ 

5k 



5k 



5E_ 

5k 



--: 5E 



(21) 



where the variation 5k is defined in terms of the functional F as in equation (fT3]) . namely, 



5k = — £ , N j ukl 



(22) 



That is, all steady states (for the evolution of any functional F) will be critical points of the 
free energy with respect to variations in free energy determined from the definition of F by 
equation (122]) in terms of ^[k]. 

The main thermodynamic postulate here is that, in principle, one can determine the evo- 
lution of the system by probing in many different directions associated with global quantities 
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F[k] (for example, the moments of a probability distribution). It is only natural to assume 
that the law for the evolution of k should be independent of the choice of which quantities 
F are used to determine it. Surprisingly, this rather general sounding assumption sets severe 
restrictions on the nature of the variation 5k. In particular, 

1. The variation 5k must be linear in 5F/5k, since the left hand side of (12T1) is also linear 
in 5F/5k. 

2. The variation 5k must transform the same way as k, as it must be dual to 5F/5k. This 
introduces the mobility fi that must be of the same type as k. 

3. The variation 5k must specify a quantity at the tangent space to the space of all possible 
k. The proper geometric way to specify this quantity is through the Lie derivative £ v 
with respect to some vector field v. 

There are only two ways to specify 5k so that it obeys these three thermodynamic and geo- 
metric constraints, when we insist that only a single new physical quantity is introduced. 
Namely, 

/ 5F\ " 

either 5k = —£ v k with v = I /jl[k] o — 1 , (23) 



5k J 

(5F \ " 
k o — I . (24) 
5k J 

Both of (1231) and (1241) are consistent with all three geometric and thermodynamics require- 
ments. However, the first possibility (1231) prevents formation of measure- valued solutions in k, 
when k is chosen to be a 1-form, a 2-form or a vector field, and we know that these solutions 
must exist. In contrast, the second possibility (1241) leads to ( 1T81) . which does admit measure- 
valued solutions for an arbitrary geometric quantity k. In the remainder of the paper we shall 
choose (124)) and investigate the corresponding evolution equation (fl8l) . The alternative choice 
(123]) would have reversed the roles of k and /a[k] in the Lie derivative. 



3 GGF vortex dynamics 

An important physical system in the context of the geometric dynamics of exact two-forms is 
the Lie-Poisson system for the vorticity of an ideal Euler fluid [22]. The developments above 
produce an interesting opportunity for the addition of dissipation to ideal fluid equations. This 
opportunity arises from noticing that the dissipative diamond flows that were just derived 
could just as well be used with any type of evolution operator, not just the Eulerian partial 
time derivative. For example, if we choose the geometric order parameter k to be the exact 
two-form uj = uj ■ dS appearing as the vorticity in Euler's equations for incompressible motion 
with fluid velocity u, then the GGF equation (TTSl) with Lagrangian time derivative may be 
introduced as a modification of Euler's vorticity equation as follows, 

d t u + £ u u = -£ (Mo sE )t uj . (25) 



Euler Dissipation 
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Euler's vorticity equation is recovered when the left hand side of this equation is set equal 
to zero. This modified geometric form of vorticity dynamics supports point vortex solutions, 
requires no additional boundary conditions, and dissipates kinetic energy for the appropriate 
choices of \x and E. Equation (125|) will be derived after making a few remarks about the 
geometry of the vorticity governed by Euler's equation. 

The Lie-Poisson bracket for vorticity is written on the dual X* ol of the Lie-algebra X vo \ of 
volume-preserving diffeomorphisms, which is isomorphic to the set of exact one-forms uo = da, 
where a is a generic one-form. In this case the Jacobi-Lie bracket between two volume- 
preserving vector fields £ x and £ 2 i n ^ 3 ma y be written as 

[£i> £2] = -curi(£i x £2) ■ 
In terms of the vector potentials for which £ x = curl and £ 2 = curl V'2 this bracket becomes 

[ £l ) £2 ] = — cur l ( cur l "01 X CUrl "02 ) • 

The vector potentials ip x and ij) 2 are defined up to a gradient of a scalar function so that 
one can always choose a gauge in which divi/> = 0. Pairing the vector field given by the Lie 
bracket with a one-form (density) a, then yields, after an integration by parts, 

("> [€1, £2]) = -(curl a, (curl Vi x curl ^ 2 ) ) = -(curl a, [[V>i, ^ 2 ]] ) 
where we have defined 

[^n V> 2 ]] := curl V>i x curl W 

The operation [[• , ■ ]] defines a Lie algebra structure on the space of vector potentials whose 
dual space may be naturally identified with exact two-forms u = curia. At this point, the 
expression for the Lie-Poisson bracket for functionals of vorticity may be introduced as 

= <«■ ttss- foil) = J u ■ ( curl ^ xcurl fc) d x ' 

where H = | J r uj ■ (—A)~ 1 ojd 3 x is the fluid's kinetic energy expressed in terms of vorticity. 

Now, vorticity dynamics is an example of geodesic motion on a (infinite-dimensional) Lie 
group [23] 

d t u) = {uj, H} = — &d* SH / Soj uj = — curl(u;xcurl (— A)" 1 ^;) = curl (curl iftx u>) = — £ CUT i^ uj 

and allows singular solutions in the form of vortex filaments. 

In order to write the GGF evolution equation (Tl8|) for u> one must compute the diamond 
operation o for the ad* action, which is defined in terms of Lie derivative by 

ad^u; = £ cxsx i^uj. (26) 

The computation of the o operation follows from its definition in equation (j3J). For any two 
velocity vector potentials (f> and if), and an exact two form uj one finds 

(0ow, V) = - (07 £ C wiii,u) = (0,curl(w x curlV)) (27) 
= ( curl <fi x uj, curl 1/? ) = ( curl (curl <f> x uj) , ijj ) . (28) 
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Consequently, up to addition of a gradient, the diamond operation is given in vector form as 

(f> o uj = curl (curl <f> x u>) . (29) 
The insertion of this expression in the bracket (120]) gives the GGF equation for u>, 



d t u) = curl ^uj x curl curl MLi[u>] x curl - — J J . (30) 

Consequently, equation fl2"5]) emerges in the equivalent forms, 

d t u = - ad^ u; + ad* (ad * mM) , w (31) 

= curl ^uj x curl ^ — + curl x curl ^ J J J . (32) 

The full dynamics for the vorticity in equation fl25l) is specified up to the choices of the mobility 
n[u}] and the energy in the dissipative bracket i?[u;]. By definition, the mobility belongs to 
the dual space of volume-preserving vector fields which is here identified with exact two-forms, 
thus one can write the mobility in terms of its vector potential as /x = curl A and rewrite the 
GGF equation (|30|) as 

d t uj = curl ( uj x curl curl A , - — 
V LL bu) 

This equation raises questions concerning the dynamics of vortex filaments with nonlo- 
cal dissipation, following the ideas in [ID] , where connections were established between the 
Marsden-Weinstein bracket [22] and the Rasetti-Regge bracket for vortex dynamics [2U 125] . 
Ideas for dissipative bracket descriptions in fluids have been introduced previously, see Bloch 
et al. [151 HE] and references therein. In particular, equation fl3"2"|) recovers equations (2.2-2.3) 
of Vallis et al. [17] when E = H and n = au) for a constant a. 



(33) 



Remark 2 The GGF equation [2by may be expressed as the Lie-derivative relation for con- 
servation of vorticity flux, 

d t (w ■ dS) = - £ u+v (w • dS) , (34) 

in which the velocities u and v may be written in terms of the commutator [ ■ , • ] of diver- 
genceless vector fields as, 

5U 5E 
u = curl - — , v = — curl curl (/4u>] x u) = curl | fJ-Ujj] , u 1 where u = curl - — , (35) 

OIjJ OUJ 

Since both u and v are divergenceless, the vorticity equation (34\ ) may also be expressed as a 
commutator of divergenceless vector fields, denoted as [•, •], 

d t u + (u + v) ■ Vw - uj ■ V(u + v) = d t v + [u + v, uj] = . (36) 

Thus, the vorticity is advected by the total velocity (u + v) and is stretched by the total velocity 
gradient. In this form one recognizes that the singular vortex filament solutions of l[3b)) will 
move with the total velocity (u + v), instead of the Biot-Savart velocity (u = curP 1 ^;) alone. 

Note that adding geometric dissipation to the vorticity equation as in equation IfSh)) does not 
require an additional boundary condition for fixed boundaries. 
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Proposition 3 The GGF equation I[2d*\) for vorticity including both inertia and dissipation 
takes the same form as the Euler vorticity equation in two dimensions, but with a modified 
stream function. 

Proof. By a standard calculation with stream functions in two dimensions, equations (1351) 
and (I34p imply the following dynamics, expressed in terms of ui := z ■ uj and fi := z ■ fi 

d t uj+ [u, ^-[/iM,^]] =0, (37) 

where ip = 5H/5uj, ip = 5E/5u and [/, g] is the symplectic bracket, given for motion in the 
(x, y) plane by the two-dimensional Jacobian determinant, 

[/, g] dx Ady = df Adg. (38) 

Equation (]37|) takes the same form as Euler's equation for vorticity, but with a modified 
stream function, now given by the sum ip — [/x, ip~\. This proves the proposition. ■ 

Remark 4 The GGF equation for vorticity in two dimensions (31) recovers equation (4-3) of 
Vallis et al. fT 7 ^ when one chooses fi = au for a constant a and E = | J ujipdxdy. However, 
for this choice of mobility, \i, point vortex solutions are excluded. 



Proposition 5 The GGF equation for vorticity in two dimensions ( 31) possesses point vortex 
solutions, with any choices of and tp for which K — tjj — \ji \uj] , is sufficiently smooth. 

Proof. Pairing equation (1371) with a stream function 77 yields 

(77, d t uj) = ( [77, /£"[£*;]], uj) where K [u] = i> - [pi [u] , $] . (39) 

Inserting the expression 

u(x, y, t) = T(t) 5(x - X(t)) S(y - Y(t)) 

into the previous equation and integrating against a smooth test function yields 

• _, • drj ■ drj dr] dK dK drj 

Tv + TX dx +TY W = T dxW- T dxW> 

where rj and K are evaluated at the point (x, y) = (X(t),Y(t)). Thus, the point vortex 
solutions for equation (1371) on the (X, Y) plane satisfy 

F^O, *.« ,40, 

whose solutions exist provided the function K is sufficiently smooth. ■ 

Remark 6 Solutions of the symplectic Hamiltonian system ( fjTj ) extend for the case of evolu- 
tion of arbitrary many point vortices for the GGF vorticity equation (31\ ) in two dimensions. 
These solutions represent a set of N vortices at positions (X^(t), Yk(t)) (k = 1, . . . , N) moving 
in the plane. Properties of the corresponding point vortex solutions of Euler's equations in the 
plane are discussed for example in JBj. 
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Remark 7 (Preservation of fluid vorticity properties in 3D) 

As a consequence of the modified vorticity equation i fffS)) in commutator form, one easily checks 
the following properties. 

• ErteVs theorem is satisfied by the vector field uj ■ V associated to vorticity. By using 
the commutator notation and the material derivative D/ Dt, one can write 



Da da , . _ _ 

:= — — h (u + v) • Va = UJ ■ Va , so that 

Dt at 

for any scalar function a(x, t). 



— , u ■ V 

Dt ' 



a = 0, (41) 



An analogue of the Kelvin's circulation theorem holds for equation f/pon expressing 

the vorticity as u) = curl u ; one writes the following dissipative form of the Euler equation 
for the velocity u 

d t u + (u + v) ■ Vu + UjVv j = -Vp , V ■ u = , (42) 
where v is given in /i35\). This equation may also be expressed as 

d t u + u ■ Vu + V (p — u ■ v J = v x curlu^ , V ■ u = , (43) 

Vortex force 

by using a vector identity. An equivalent alternative is the Lie derivative form, 

(d t + £ u+v )(u-dx) = -dp. (44) 
Hence, we find that a modified circulation theorem is satisfied, 



d 

Tit 



u ■ dx = (45) 

C(u+v) 



for a loop C(u + v) moving with the "total" velocity u + v. That is, two velocities appear 
in the modified circulation theorem. One is the "transport velocity" u + v and the other 
is the "transported velocity" u. 

• From equations ( [73] ) and (3$ ) one checks that 

(d t + £ u+v ) (w ■ dS A u ■ dx) = - w ■ dS A dp = -div (p w) d 3 x (46) 

so that the helicity of the vorticity uj is conserved 

*JJL u - u * x=o (47) 

We may summarize these remarks as follows: 

All of these classical geometric results for ideal incompressible fluid mechanics 
follow for the modified Euler equation. These results all persist (including preser- 
vation of helicity) when transport velocity is replaced as (u + v) — > v. 
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Remark 8 (Relation of the GGF vorticity equations to Craik-Leibovich theory) 

The Craik-Leibovich (CL) equations describe the dynamics of the Eulerian mean fluid ve- 
locity u depending on time t and spatial position x in three dimensions, when the fluid motion 
is driven by rapidly oscillating surface waves due to the wind. These circumstances may gen- 
erate Langmuir circulations - sets of vortices with axes nearly parallel to the wind direction 
which sometimes occur in the upper layers of lakes and oceans. 

In the CL theory, the rapidly oscillating waves at the surface are assumed to be unaffected 
by the more slowly changing currents below. The effect of the waves on the Eulerian mean flow 
is parameterized in the CL theory by introducing into the Navier-Stokes equations a "vortex 
force," expressed in terms of a prescribed Stokes drift velocity, ug(x.,t). 

The CL equations are given by, 
du 

— + (u • V) u + = u s x curl u + uAu , V ■ u = , (48) 

at v v ' 

Vortex force 

where w is a pressure enforcing incompressibility and v is viscosity, ignored hereafter. Ones 
sees that the GGF vorticity equation j^M ) is in the same form as the Craik-Leibovich (CL) 
equation §JjM . However, the velocities in the vortex force in the two cases differ. In both cases, 
the vortex force may be interpreted as a noninertial force arising from having transformed into 
a prescribed moving frame. For additional references and discussions of the properties of the 
CL equations, see also [27^ . 

This completes our investigation of GGF vortex dynamics. An obvious extension would be 
to consider GGF vortex patches in two dimensions. For example, we leave to another study the 
investigation of the selective decay hypothesis of Matthaeus and Montgomery [28] for the 
relaxation to vortex equilibrium states. Instead of pursuing such GGF vorticity considerations 
further, we shall next turn our attention to the geometric gradient-flow (GGF) equation for 
an arbitrary geometric quantity. In particular, we shall perform explicit computations of GGF 
equations for scalars, one-forms and two-forms. In each case, the GGF evolution equations 
allow singular (measure-valued) solutions. Thus, the connecting theme of the discussions 
below with the Euler equations is the presence of singular solutions in GGF equations, akin to 
point vortices, vortex lines or vortex sheets for the Euler equations, but with geometric order 
parameters that transform differently from vorticity under smooth invertible maps. 



4 Singular GGF solutions for scalars 

We have seen in section 12.21 how GGF equations (I18p always allow for singular solutions of 
the form 

N 

k(x, t) = 22 P<*( s > t)S(x- q a (s, t)) ds . (49) 

a=0 J 

To derive the equations for p(s,t) and q(s,t), we substitute the solution ansatz (j4"9l) into the 
GGF equation ffTB"]) and integrate it against a smooth test function <f> dual to k, thereby obtain- 
ing (JHD- This section illustrates one particular example. Two classes of geometric quantities 
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admitting singular solutions of (|18[) are known [12]. One class includes, for example, scalars, 
1-forms and 2-forms, and gives characteristic equations for which the characteristic velocity is 
a nonlocal vector function. A second class, which includes densities (three-forms on M 3 ), is a 
nonlinear nonlocal gradient flow equation (flTT) . Since our main interest lies with the first case, 
we shall concentrate on the nonlocal characteristic equations. These characteristic equations 
have interesting mathematical and physical properties and many potential applications. 

The fundamental example is an active scalar, for which k = f is a function. The evolution 
of a scalar by equation ffT8l) obeys 

^/ = -%/]<>f)»/ = -(^ V M/]) fl - V/ - (50) 
Equation (150]) can be rewritten in characteristic form as 

{5E \» 

df/dt = along dx/dt = (jj^M) ■ ( 51 ) 

The characteristic speeds of this equation are nonlocal when SE/Sf and \i are chosen to depend 
on the average value, f. It is interesting that such problems arise commonly in the theory of 
quasi-geostrophic convection and may lead to the development of singularities in finite time 

|2Sl EDI EI]. 

Explicit equations for the evolution of strengths p a and coordinates q a for a sum of 8- 
functions in (1491) may be derived using (jSJ) when fi[f] = H * f = f. The singular solution 
parameters satisfy 



dp a (t,s) (8E \t 

— — — = p a (t, s) div \— V/4/j 1 



(52) 

=«<,(*>«) 



dqjf, s) (5E \t 

Pa(t,s) — = Pa{t,s)^—Vp J [f]J 



(53) 

=</a(M) 



for a = 1, 2, . . . , N. For the choice fj, [/] = /, a solution containing a single 5-function satisfies 
p = — Ap 3 , so an initial condition p(0) = po, evolves according to l/p(t) 2 = 1/pq + AaH. The 
comparison of 1/p 2 from numerics with this theoretical prediction is shown in Fig. [TJ 



5 Active one-forms and two-forms 

In more generality, one may develop a nonlocal characteristic equation for the evolution of any 
geometric quantity. We write the explicit form of the equations for the examples of differential 
1-forms A • dx and 2-forms B • dS in three-dimensional space. For this, we begin by computing 
the the Lie derivative and the diamond operation for these cases. In Euclidean coordinates, 
the Lie derivatives for these two choices of k are: 

- £ v (A • dx) = - ((v • V)A + AjVV) • dx 

= (v x curl A — V(v ■ A)) • dx , 

-£ v (B-dS) = - d(v J (B • dS)) - v J d(B ■ dS) 

= -d((v xB) ■ dx) - v _l(dwBd 3 x) 

= (curl (v x B) — v div B) • dS 
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Figure 1: Evolution of the ^-function strength l/p(t) 2 versus time (circles) from [12]. The 
theoretical prediction 1 /py+At is shown as a solid line obtained without any fitting parameters. 



Both of these expressions are familiar from fluid dynamics, particularly magnetohydrodynam- 
ics (MHD). 

From these formulas for Lie derivative in vector form and the definition of diamond in 
equation (j3J), we compute explicit expressions for the diamond operation with 1-forms and 
2-forms, 

FA1 6E \ / 5E , f A 1 r « l j . SE 

^ * °6A ,U ) = \6A X ^ * ~ ™ * 6A ' U 

SE \ I , 5E 5E , 

M[B]o — , U/ ) = (fj,[B\ x curl— - — div/z[B], u 

for any vector field u. 

Remark 9 These formulas result from computing the diamond operation using the standard 
action of Lie derivatives on differential forms, instead of the action used earlier to calculate 
equation l29\) which employed vector potentials for divergenceless vector fields. 

The explicit forms of these equations are unfamiliar and their solutions will be investigated 
elsewhere in the context of their potential applications. Here, we discuss the cases of exact 



18 



one- and two-forms. In these cases curl/x[A] = = div//[B], so the equations simplify 
considerably. In vector notation one now has 



dA 
~dt 



Vn/i[A]div^y. Aj (54) 



<9B / / 5E S 

curl fi[B] x curl — 1 X B I (55) 



dt \ V 1 1 5B 

Now, upon defining the vector field 



v x 



( r , <^Y 
ll\k\ o — with k, = A, B 



the following result holds. 



Proposition 10 The geometric gradient-flow equations (54\ ) and 153]) for closed one-forms 
A two-forms B have singular solutions of the form (f^ffi ), where 

q a (t,s) = v(a;)| 3;=qa 

Pa(*,s) = p a (t,s) (V-VCX))!^ - Vv^U^- Po (t, S ) (56) 

for closed one-forms A and 



q a (t,s) = v(x) 



x =<i a 



pa(t, S) = p a (t, S) ■ Vv(cc) 

for closed two-forms B. 

Proof. Consider equation (j54p for A written as 

fi div — ] 

5AJ 

Pairing this equation with a smooth vector field 0, substituting the singular solution ansatz 
(f49|) . integrating by parts where necessary and matching all terms in (p on the two sides yields 
equations (|56l) for the q's and p's. 

The result for closed 2-forms is proven by noticing that 

fi x curl — 

then following the same steps as for the case of exact 1-forms. ■ 

When considering the GGF equations ffl8|) for differential forms that are not exact, singular 
solutions also exist, satisfying more complicated relations, depending on how the diamond op- 
erator is expressed in each particular case. One may see this, by following the same procedure. 
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For exact one- and two-forms, the vector equations above can be reduced to nonlocal 
nonlinear scalar characteristic equations of the form (IBUj) for the potentials:. Note that in IR 3 
(which is of interest to us here) every closed form is exact since curl A = gives A = for 
some scalar ip and div B = necessitates B = curl C for some vector C. The characteristic 
equations for the potentials are derived in the following 

Proposition 11 The vector equations (54\ ) and |53]) for exact 1-forms A = and exact 
2-forms B = curl {m z) are equivalent to scalar GGF equations of the type [5fy) , in terms of 
the potentials ip and m . Specifically, one finds 

£-(«fv*i)'.V,. (57, 



\8ip 

V$[tf] ) ■ W, (5* 



and 

dm _ /5£_ Tr _\» 
~dt ~ \5m 

where one defines jtt[A] := Vi9[^] and /x[B] := curl (&[m] z). 
Proof. Inserting the expression A = Vip in eq. ( 1541) yields 

dip /fait SE\i _ , 
/ 5E\I 

= [vm —) • 

with nonlocal 5E/b~ip and /J,[ip]. 

Similarly, the evolution of 2-form fluxes B ■ dS = B x dy A dz + B y dz A dx + B z dx A dy also 
simplifies, when B = Vm x z where \l/ only depends on two spatial coordinates (x, y). Then, 

JE JE 
curl — = z— . 
5B dm 

and 

, SE ^ „ SE SE _ 

uB x curl — = (V$ x z) x z — = - — V$ . 

Equation ( 1501) may be written for the stream function \1/ (removing the curl from both sides 

of J55D) 

*m = ~ (m v *) x b < 59 > 

Then, simplification of two cross products leads to 

_=(_V.).W. (60) 

Hence, choosing 5E/5m and $ to depend on the average value m again yields a nonlocal 
characteristic equation. ■ 
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Remark 12 Equations (51) and l[58\) do allow singular 5 -like solutions of the form fifty for ip 
and These solutions, however, lead to 8' -like singularities in the forms A and B. One may 
understand this point by deriving the expressions for ip and \1/ corresponding to the clumpon 
solutions of the form ( [7^ ] for A and B. 

For example, taking the divergence of an exact one-form A = Vip yields V ■ A = Aip. 
Upon using the Green's function of the Laplace operator G(x, y) = — \x — y \~ , an expression 
forip emerges in terms of A: 

tp(x,t) = - J V x >G(x,x') -A(x',t) dx. 

Inserting the singular solution ( [7^ for A then yields 

il>{x, t) = -J2 [ds Pi{s, t) • V Qi G(x, Q^s, t)). 

However, this singular solution for the potential is not in the same form as ffllfy , since the 
singularities for ip do not manifest themselves as 5-functions. 

A similar procedure applies to the case of exact two-forms B(sc, y) = curl(^f(x,y) z), so 
that curlB = A z) . We have 

*(x) = z • [ dsPi(s,t) x V Qi G(x,Qi(s,t)), 

where Q is in the plane (x, y) . 

Thus, the equations (57| ) and / f5^) for ip and allow for two species of singular solutions. 
One of them takes the form fijjty , while the other corresponds to a 5-like solution of the same 
form ( [^PP for A and B. 

A deeper explanation of this fact can be given in a general context as follows. Consider 
the advection equation for an exact form k = dX, with potential A 

(dt + £ u ) d\ = 0. 

At this point, one remembers that the exterior differential commutes with the Lie derivative so 
that the equation for the potential A is again an advection equation with the same characteristic 
velocity 

(d t + £ u ) A = 

At this point, one obtains singular 5-like solutions of the form ftjlty) for both k and A (provided 
the characteristic velocity u is sufficiently smooth). 



6 Three more examples 

The developments discussed above produce an interesting opportunity for the addition of dis- 
sipation to various other continuum equations. Following the introduction of the dissipative 
Euler equation above, one could extend the dissipative diamond flows with any type of evolu- 
tion operator. This section sketches how one might develop this idea further, by illustrating 
its application in three more physically relevant examples. 
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Dissipative EPDiff Consider adding geometric dissipation to the Euler-Poincare equation 
on the diffeomorphisms (EPDiff) [5] for the evolution of a one-form density m defined by 

m = m • dx (g> d 3 x . (61) 

This addition gives the dissipative EP equation, 

d t m + &dg H/Sm m= - £ f u m = -&d* m . (62) 

' yn[m]o8E/8m) U[m] o6E/Sm J 

When H[m] is the Hamiltonian for the Lie-Poisson theory corresponding to EPDiff, then the 
vector field SH/Sm = u is the characteristic velocity for the Euler-Poincare equation. For 
a one-form density m, the diamond operation is given by ad*, which is equivalent to Lie 
derivative. That is, 

fi[m] o 5E/Sm = &d* SE/Sm fi[m] = £ 6E/6m ^[m} . (63) 

The further choice /i[m] = am for a positive constant a recovers equation (6.10) of Bloch 
et al. [15] . When in addition, /z[m] = K * m for a smoothing kernel K, then equation (!62l 
supports singular solutions of the type discussed in Holm and Marsden [3]. 

Peakon dynamics for the GGF modified Camassa-Holm equation In one dimension, 
the GGF version of the EPDiff equation (162]) reduces to, 

d t m + (u + v)m x + 2m(u + v) x = , (64) 

where u = 5H/5m for a specified Hamiltonian H[m}. The other velocity v is given in one 
dimension by 

v = ( &d* 5 E/5m ) " = ^ d x fj,[m) + 2/i[m] ^ — , (65) 

for arbitrary (smooth) choices of fi[m] and i?[m]. Now consider the singular solution form for 
m given by a sum of N delta functions, 

N 

C, t) = 

i=l 



:,t)=X)Pi(*)^ -?<(<)), (66) 
i=i 

and take quadratic functionals if [m] = 1/2 (m, G * m) and i?[m] = 1/2 (m, VF* m) so that 



and 



N 

u(x) = G * m = pj G(x — qj) 

3=1 



v(x) = W* m d x K * m + 2 K * m d x W* m 

N . 

= p i Pk [ K ( X ~ q ^ d * W ( x ~ Ik) + 2 W(x - q k ) d x K(x - qj) 

N 

j,k=i 
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where we have defined 1Z for compactness of notation. Substituting the above expressions 
into the GGF EPDiff equation (I6"4l and integrating against a smooth test function yields the 
following relations for time derivatives of Pi(t) and qi(t): 



u{x) + v{x t 

N 



x=qi(t) 



(67) 



N 



Pi 



Pj G ^ ~ 3j) + X Pi P k ^(ft ~ 3j> ft 
j=l J,fc=l 

— Pi iu'(x) + v'{x) 



qkj 



\x=qi(t) 

N N 

Pi Y Pi -Qj) ~ Pi X Pfc 9 <z^(ft - 1i, Qi 



(6* 



<?fcj 



i=i j,fc=i 

The choices of H, E and /x as functionals of m determine the ensuing dynamics of the singular 
solutions. In particular, in the case when the velocity u is given by 



u\m\ 



5H 

— - = (1 - a a, 
dm 



*m = - I e 
2 



J e ' x x '^ a m(x')dx' , for H = ^J mu[m]dx. (69) 



Equation (1641) is a GGF version of the integrable Camassa-Holm equation with peaked soliton 
solutions pE]. Nonlinear interactions of N traveling waves of this system may be investigated 
by following the approach of Fringer and Holm 



Dissipative Vlasov dynamics One may also extend the diamond dissipation framework to 
systems such as the Vlasov equation in the symplectic framework of coordinates and momenta 
as independent variables. This extension requires the introduction of the Vlasov Lie-Poisson 
bracket, defined for phase space densities f(q,p) dq A dp on T*M N as 



{F,H}(f) = J! f 



SF 5H 
Jf' Jf 



dq A dp . 



(70) 



QP 



Here the Lie bracket [• , - \ qp is now given by the canonical Poisson bracket on Hamiltonian 
functions. At this point a new definition of the diamond operator is required. This is found 
by the well known identification of Hamiltonian vector fields and their generating functions 
that gives the relation 

-£x h f = \f,h] v (71) 

for any phase space density f(q,p) and any function h(q,p). This relation identifies the 
symplectic Lie algebra action on the Vlasov distribution and the new kind of symplectic 
diamond (which we denote by *) which may be computed by applying the general definition 

as 

ig*f,h) = {g,-£ Xh f) = ([g,f] qp ,h) (72) 

for any two functions g and h. Extending the previous discussions, one can write the following 
form of GGF dissipative Vlasov equation, 



dt 













\f - 


qp 


f, 




qp _ 



(73) 



23 



where, in general, the functionals H and E are independent. This equation has the same form 
as the equations for a dissipative class of Vlasov plasmas in astrophysics, proposed by Kandrup 
[33] to model gravitational radiation reaction. Kandrup's formulation for an azimuthally 
symmetric particle distribution is recovered by choosing a linear phase space mobility u = af 
with positive constant a and taking E to be J z [f] the total azimuthal angular momentum for 
the Vlasov distribution /. More generally, if one chooses u[f] = af and E to be the Vlasov 
Hamiltonian Hf, the dissipative Vlasov equation (1731) assumes the double bracket form, 



df 
dt 



+ 



/ 



5H 



f 



Sf 



a 



<!!> 



f 



f 



5H 



f 



Sf 



(74) 



<jp 



This is also the Vlasov-Poisson equation in Bloch et al. [15]. However, in contrast to the 
choices in [151 HHJ [2QJ [33], the GGF form of the Vlasov equation (1731) allows more general 
mobilities such as u[f] = K * f (which denotes convolution of / with a smoothing kernel 
K) . The GGF choice has the advantage of recovering the one-particle solution as its singular 
solution. 



Dissipative semidirect product dynamics The equations derived above consolidate the 
idea that any continuum equation in characteristic form, 

(d t + £ u ) « = , 

may be modified to include dissipation via the substitution u — > u + v, in which v is the 
dissipative velocity term expressed in equation (1T81) . This idea may also be extended to the 
semidirect product framework presented in [M], in order to include compressible fluid flows 
and plasma fluid models such as MHD. Applications of the semidirect product framework 
for continuum mechanics are beyond the scope of the present work. However, it would be 
interesting to pursue these applications in future investigations. An immediate application of 
the semidirect product framework would be an investigation of ideas of selective decay in 
the approach to topological equilibria for example in MHD, as first suggested by Taylor 
[55] based on work of Woltjer [SB] and later elaborated by Moffat [H7J and others. A possible 
counterpoint would be to treat the additional helicity- conserving geometric force as a means 
of driving a magnetic dynamo, rather than relaxing to an equilibrium. 



7 Conclusions and Comments 

This paper has provided a general method of constructing dissipative evolutionary equations 
in the form (1TT1) for a variety of different types of geometric order parameters n. The method 
produces a plethora of fascinating singular solutions for these evolutionary GGF equations. 
Each GGF equation is expressed as a Lie derivative (TTT1) . Thus, these are all characteristic 
equations in a certain geometric sense. However, the characteristic velocities in these equations 
may be nonlocal. That is, the characteristic velocities may depend on the solution in the 
entire domain. The equations may possess either or both of the following structures: (i) a 
conservative Lie-Poisson Hamiltonian structure; (ii) a dissipative Riemannian metric structure. 
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The Hamiltonian evolution of a continuum material by the Poisson structure is similar to 
Lagrangian fluid dynamics, while the dissipative evolution by the Riemannian metric structure 
is a geometric version of gradient flow similar to Darcy's Law. The two types of evolution are 
combined by simply adding the characteristic velocities in their Lie derivatives. 

Similar types of equations were discussed by Bloch et al. [15] who studied the effects on 
the stability of equilibrium solutions of continuum Lie-Poisson Hamiltonian systems of adding 
a type of geometric dissipation that preserves the coadjoint orbits of the Hamiltonian systems. 
Such equations have the form 



for two bracket operations, one antisymmetric and Poisson ({F, H}) and the other symmetric 
and Leibnitz ({{F,H}}). In contrast, the GGF theory applies to a more general family of 
such equations with two independent types of energy H and E and two types of brackets, one 
antisymmetric and Lie-Poisson, the other derived in (|2"U1) from thermodynamic principles; so 
that systems of the GGF form 



need not preserve coadjoint orbits, but can have the additional feature that it admit singular 
solution behavior. The GGF theory has been shown to apply in a number of continuum 
flows with geometric order parameters, in examples ranging from gradient flows governed by 
a nonlocal version of Darcy's Law, to dissipative modifications of Euler fluid flows and plasma 
dynamics. The various types of singular solutions include point vortices, vortex filaments and 
sheets, solitons and single particle solutions for Vlasov dynamics. These singular solutions 
apparently also apply to a recent generalization of Landau's two fluid model of superfluid 



Each of these equations admits singular solutions which lie on invariant solution manifolds. 
In some cases, the singular solutions emerge from smooth confined initial conditions. In other 
cases, such emergent behavior does not occur. It remains an open question to determine 
whether the singular solutions of a given geometric type will emerge from smooth initial 
conditions. After they are created, the singular solutions evolve with their own dynamics. 
Investigations of the interactions of these singular solutions and the types of motions available 
to them will be discussed elsewhere. The present paper has derived the dynamical equations 
for these singular solutions in various cases. However, investigations are left to the future to 
characterize their formation, stability and interactions, the integrability of their interaction 
dynamics and their fundamental mathematical nature. 
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